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ABSTRACT 

This paper reports on a viscoelastic model for deformation and stress associ- 
ated with earthquakes. The model consists of a rectangular dislocation (strike-slip 
fault) in a viscoelastic layer (lithosphere) lying over a viscoelastic half-space (as- 
thenosphcrc). The first part of the paper contains an analysis of the time-dependent 
surface stresses. The model predicts that near the fault a significant fraction of the 
stress that was reduced during the earthquake may be recovered by viscoelastic 
softening of the lithosphere. By contrast, the strain shows very little change near 
the fault. The model also predicts that the stress changes associated with asthen- 
ospheric flow extend over a broader region than those associated with lithospheric 
relaxation even though the peak value is less. The second part of the paper studies 
the dependence of the displacements, strains, and stresses on fault parameters. 

Peak values of strain and stress drop increase with increasing fault height and de- 
crease with fault depth. Under many circumstances postseismic strains and stresses 
show an increase with decreasing depth to the lithosphcrc-asthenospherc boundary. 
Values of the strain and stress at distant points from the fault increase with fault 
area but are relatively insensitive to fault depth. 


P0STSB1SMIC VISCOELASTIC DEFORMATION AND STRESS 
Part 2: Stress Theory and Computation; Dependence of 
Displacement, Strain, and Stress on Fault Parameters 


I. INTRODUCTION 

This paper is devoted first to a consideration of the postscismic stresses associated with earth- 
quakes using a viscoelastic model of the relaxation processes. It also considers the dependence of 
the displacements, strains, and stresses on fault parameters. The mathematical model consists of a 
vertical rectangular dislocation located in a viscoelastic layer which lies over a viscoelastic half- 
space with different rheological properties. Figure 1 shows the relevant features of this model 
which was used in previous work (Cohen, 1979b; herein referred to as Part 1. Sec also Cohen, 1979a) 
to calculate postscismic displacements and strains due to earthquakes along strike-slip faults. The 
upper layer, of thickness II, is meant to represent the lithosphere. Its dcviatoric rheology is 
modeled by a standard viscoelastic solid; its ditntatiomil rheology will generally be taken to be 
clastic (model D1LE), although we have also considered the standard viscoelastic solid case. The 
half-space represents the asthcnosphcrc for which only the dcviatoric response needs to be con- 
sidered. This we have taken to be a Maxwell substance with a response time much longer than 
that of the lithosphere. The dislocation which models the fault has length, 2L, is located at a 
mean depth, 5, has a height, Ah. The size of the dislocation at the time of the earthquake is U 0 . 
Additional details about this model were discussed in Part 1 where we also presented the computed 
displacements and strains. These strains arc used in the present paper to compute coseismic and 
postscismic stresses taking into account the assumed viscoelastic properties of the earth. The 
analysis and numerical results will indicate that the patterns of postscismic stress and strain may 
be considerably different. In particular the postscismic shear stress recovery near the fault follow- 
ing an earthquake may be large even when tile strain change is small. Thus postseisinic deforma- 
tions may be due to both viscoelastic changes in rigidity and postseismic fault motion caused by 
stress recovery. 
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II. STRESS CONSIDERATIONS 


The evaluation of the coseismie and postseismic stresses starts with a consideration of the 
stress, o, - strain, e, equation for a homogeneous isotropic elastic body, viz., 

ff ij = 2/*i ey + (k, - | ^ A fiy (1) 

wiiere jaj and kj arc the shear and bulk modulus, A is the dilatation (A = e n + e 22 + e 33 ) and 
5jj = 1 if i — j or 5jj = 0 if i =£ j. To determine the surface stress for viscoelastic material \vc adopt 
the following procedure. First attention is focused on three times. The first of these is t 0 , the 
time immediately after the earthquake; the second tg is a time long compared to the lithosphere 
viscoelastic relaxation time and short compared to the asthenosphere relaxation time; the tiiird t a 
is a time long compared to the asthenosphere relaxation time. The viscoelastic strains at cadi of 
these were calculated in Part I where we argued that the viscoelastic displacement equations can 
be obtained from a knowledge of the clastic displacement equations by replacing the elastic moduli 
by moduli appropriate to the viscoelastic case. Similarly if the viscoelastic strains arc known then 
the correspondence principle (Flugge, 1967) and the viscoelastic constitutive equations indicate 
that the stresses may be computed from Equation (1) provided n 1 and k 2 arc replaced by the 
appropriate effective moduli. Specifically the constitutive equation for the deviatoric behavior of 
the lithosphere is 


o + 


u t n 

Mn + Mb M a + Mb M a + Mb 


( 2 ) 


At t 0 there is a sudden change in stress and strain so terms in a and e are large compared to those 
in a and e. Tims 

Ao(t 0 ) = ju a Ae(t 0 ) (3) 

Tims the effective clastic shear modulus of the lithosphere at time t 0 is, from the correspondence 


principle, 


, _ Mn 
M] (to) “ ~ 


( 4 ) 
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For times tj» and t ;l the terms in a and e dominate the terms involving the time derivatives so 


a(tj> or t :l ) = 


P u 

1 + j3 


e(t e or t a ) 


where ft = In this case 


M tfi°rt a ) = — 


Pa 


2 1+0 


The dilatior.nl stress-strain law of tlie litliosplicre is 


a = ke 


So from tile correspondence principle (and at all times) 


( 5 ) 


( 6 ) 


(7) 


k , ■!- 


( 8 ) 


Witli these results we find 


auCto) 1 / k \ 

+-(— - 1 A(t 0 )S,j 

M a 3\n a / 


gjjCtfi) _ ejjCtju) 1 / k 1 

— ' i ~+7 + 1 w" r +'j3. 
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gjj^a) = gjjftn) + J_ A 
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1 + P 3 \ /u., 1+0 


|A(t c )5ij 


|A(t a )5ij 


(9a) 


(9b) 


(9c) 


k|(t 0 ) 2 k 5 

In our numerical computations we have taken = — corresponding to a Poisson’s 

Pi ( f o ) 3 Pu ^ 

ratio of 0.25 and have assumed 0 = 3/2. The evaluation of the surface shear stress, a 12 , is ob- 
tained directly from Equation (1) using the shear strains computed in Part 1 with the assumption 
that the instantaneous shear moduli of the lithosphere and asthenospherc are equal. The evalua- 
tion of the surface normal stresses, o n and a 22 , requires knowledge of the dilatation A = + 

e 22 + c 33 - The horizontal strains, e u and e 22 , were evaluated in Part 1 under the condition that 
at the earth’s surface o 33 - 0. Inserting this condition in Equation (9) permits evaluation of e 33 
(and hence A, and o^) in terms of e n and e 22 . Explicitly 
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(10a) 
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til (to) + e 22 (*())] 
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k 
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[ e ll (tfi or t a ) + e 22 (tfi or t a )] 


(10b) 


III. STRESS COMPUTATIONS 

Tlie earthquake related stresses have been evaluated at the tiirce times t 0> tjj, and t a . The 
results shown in Figures 2-4 are for the same conditions used in computing the displacements and 
strains in Part I namely L = H = 4Ah = 86. The shear stresses a 12 (t 0 ), a I2 (t fi ), and o I2 (t a ) are 
proportional to the corresponding shear strains with the proportionality constant being n a at t 0 

jJL a 

and at to and t„. The contribution to the stress due to lithospheric relaxation can be 

1 + j3 

analyzed by defining a shear modulus difference by 


a - _ 0^a 

~ ” — Pq - — 


1 + 0 


1 + 0 


01 ) 


Then 


~ ^(Ifi) ~ Ae JZ + A/ie I2 (to) (12) 

1 + p 

where 

A<=i 2 - e^Ctg) ~ e i2^o) (12) 

Tlie change in the shear stress due to lithospheric relaxation is composed of two components. One 
contains the change in strain, tlie other the change in rigidity. Of course these components are coupled 
to one another since the strain changes are also determined by the changes in lithospheric rigidity. 
It is interesting to compare the stress changes with tlie corresponding strain changes. Such a 
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comparison is shown as a function of distance from the fault along the y axis in Figure S. At 

many points, particularly near the fault, Ae l2 is small and Ao, 2 is dominated by Afi. At the 
Aei 2 Ao n 

fault, for example, — — 10* 2 but — ~0.6. Thus there is a significant postseismic 

e i2(*o) cri 2 (t 0 ) 

stress recovery even thougii the strain change is small. Clocely related arguments concerning stress 
recovery following an earthquake have been advanced by Dieterich (1972), Cohen (1978), and 
Yamashita (1979). As for the changes occurring as the asthenosphere relaxes, they are proportional 
to the strain changes since 

[e ^ ' (14) 

i T p 

When these results arc combined with those of Part 1, the following picture of postseismic phe- 
nomena emerges. In the short term and near the fault the direct postseismic deformations due 
to rigidity changes arc small. Since, however, the stress recovery may be significant, postseismic 
deformation may be due to ascismic afterslip on the fault. The mechanisms for such afterslip may 
include stress recovery along with time dependent friction (Dietrich, 1972), space dependent fric- 
tion (Cohen, 1978), fatigue failure (Scholz, 1972), etc. Faraway from the fault the effects of 
afterslip are small unless the source is very deep. Thus on a relative basis distant postseismic de- 
formations seem to be more attributable to the softening of the lithosphere rigidity. As for the 
long term postseismic deformations the model attributes them to the flow of the asthenosphere. 

The stresses extend over a broader zone than they do in tiie case of the lithosphere relaxation. 

The relative senses of the strains and stresses due to the instantaneous response, the lithospheric 
relaxation, and the asthcnospheric flow may differ from one another and they vary with the source 
parameters and observation point. 

IV. DISPLACEMENT, STRAIN, AND STRESS VARIATIONS WITH FAULT PARAMETERS 
The calculated values of the viscoelastic displacements, strains, and stresses arc dependent on 
the choice of model parameters that arc introduced into the numerical evaluations. These param- 
eters can be grouped into two groups. The first group is comprised of the material parameters, 
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namely tlic elastic moduli and viscosities, The second group is composed of geometric factors 
including the fault dislocation, U 0 , the lithospheric thickness, H, the coordinates of the observa- 
tion point under consideration, and the parameters defining tile fault: the length, 2L, the depth, 

5, and tire height Ah. In this section sve examine some aspects of the dependence of the deforma- 
tions and stresses on the geometric factors. To restrict the discussion to a reasonable length it is 
convenient to examine how the results vary with distance from the fault along the y axis. The 
non-vanishirtg quantities along this axis arc the parallel displacements, u(y,t) the shear strains, 
e J2 (y ,t), and the shear stresses, a, 2 (y,t). Specifically, we consider the initial values: u(t 0 ), ej 2 (t 0 ), 
a ]2 (t 0 ), the differences between quantities at times t g and t Q : u(t g ) - u(t 0 ), e 12 (t g ) - e 12 (t 0 ), “ 

<7 J2 (t 0 ), and the differences between quantities at times, t a and t g : u(t a ) - u(tg), Cp(t a ) - e 12 (t g ), 
a i2( t a) - cr 12 (tj>), As before, it is convenient to normalize the results to a unit value of fault slip 
leaving as independent parameters, 5, All, II and L. The following paragraphs summarize general 
tendencies we have noted for the dependence of the peak values and distant values of deformation 
and stress on the fault parameters. Our remarks are based on calculations with parameters in the 
ranges 5<5<l5,10<Ah< 30,20 < II < 100, 50 < L<250 where all numbers arc in kilometers. 

It should be noted tiiat in some cases there are exceptions to the general tendencies we discuss. 

For ruptures tiiat penetrate the surface, the peak values (in an absolute value sense) of the follow- 
ing quantities occur at the fault: u(t 0 ), e 12 (t 0 ), a n (t Q ), ^(tg) - e 12 (t 0 ), <x 12 (t g ) - a 12 (t 0 ), e 12 (t a ) - 
G n a 12 ~°12 (tfi)- The peak values of u(t g ) - u(t 0 ) and u(t a ) - u(t g ) occur at some distance 

from tiie fault. For ruptures that do not penetrate the surface some caution must be exercised. For 
example, in tills case, cr j2 (t 0 ) is positive, indicating a stress rise, at the surface point above the fault. 

We take as a maximum value of the surface stress drop the most negative value of a J2 (t 0 ) which 
occurs at some value y>0. Similar remarks apply to some of the other quantities we discuss 
below. Variations ill 6, Ah, L, and H have the following general effects: 
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Depth, 5 


The quantities u, e 12 and aj 2 at time t 0 and the differences between these quantities at times 
tg and t 0 decrease with increasing 5 because deeper sources in the lithosphere produce smaller 
peak surface effects. The differences between these quantities at times t., and tg increase with 6 
due to the reduced distance to the asthcnnspherc with the attendant greater lithosphcrc- 
asthenosphere interaction. 

Fault Height, Ah 

All the peak displacements, suains, and stresses increase with Ah due to the increased source 

size. 

Length, L 

For ruptures that break the surface, the peak displacement occurs at the fault and is a boundary 
condition independent of L. For buried faults, the peak coscismic displacement tends to decrease 
somewhat with length. In addition the coseismic strain and stress cither decrease with or are independ- 
ent of fault length. The differences between the deformation and stress quantities at times t a and tg in- 
crease with L, presumably due to the greater lithosphere-asthenosphere coupling with the larger fault 
length. The quantity u(tg)-u(t 0 ) also increases with increasing L. By contrast c J2 (tg)~e 12 (t 0 ) de- 
creases with increasing L. So docs o 12 (tg) — o- 12 (t 0 ) but the dependence is weak in this latter case. 

Lithospheric Thickness, M 

Tile values of displacement, strain, and stress at time t 0 are independent of H due to the 
assumption p,(t 0 ) ~ MaCto)* This assumption leads to the condition that the interface between 
the lithosphere and asfhenosphcre is not sensed at the time of the earthquake and the earth re- 
sponds as an elastic half-space. By contrast, peak values of all differences between displacements, 
strains, and stresses at times t a and tg decrease with increasing H due to the increased distance 
between the fault and the lithosphere-asthenosphere boundary. The quantities u(tg) - u(t 0 ) and 
~ < 42 ^ 0 ) a f so decrease with increasing H. By contrast p 12 (tg) - a 12 (t 0 ) is affected very 
little by changes in H. 
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In a manner mimicing the preceding discussion we can discuss the dependence of the dis- 
placements, strains, and stresses of the favlt parameters at distant points from the fault. To do 
this, we must first note that at distant points (y > L, H) there are several phase changes in e 12 
and a 12 compared to their peak values. Specifically £ 12 On) and o J2 (t 0 ) are positive and Oj 2 0fi) “ 
(7 J2 Oo) is generally negative. In addition to our general comments beiuw, we note two very simple 
relationships that follow from the displacement and strain equations for y large: 


u(to) = 


e 120o) 


2LAh 


47ry 2 47ry 2 

2LAh A 


47ry 3 47ry 3 


(15) 

(16) 


where A is the area of the fault. 


Depth, 6 

To a fair approximation n, e l2 , and a l2 arc independent of depth unless 5 approaches H. 
Since at large distances from the fault the exact depth of the source matters little so long as the 
depth is shallow. The most sensitive dependence of the distant deformation variables on 5 is that 
of U(t a ) - u(f fi ). 

Fault Height, Ah 

Increasing the fault height increases u, e 12 , a J2 and the corresponding time differences. 

Fault Length, L 

Increasing the fault length increases u, e 12 - o J2 and the corresponding time differences. 
Lithospheric Thickness, I-i 

Based on the previously discussed assumption that juj (t 0 ) = ju 2 (t 0 ) the lithospheric thickness 
does not effect u(t Q ), e l2 (t 0 ), or o 12 (t 0 ). The quantities u(t B ) - u(t 0 ), e 12 (t fi ) - e J2 (t 0 ), u(t a ) - 
u(tg), e J2 (t a ) " e i 2 ( l g)> and ff i 2 ( t a) “ Uj 2 (tje) generally decrease with increasing H as do the 
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corresponding peak values. The quantity a, 2 (tg) - o, 2 (to) e ‘^ ,cr decreases with increasing H 
or shows little change depending on whether the change in n, or in e 12 is controlling the change 
in Oj2« 

V. SUMMARY 

This paper has considered a model of the surface stresses associated with viscoelastic relaxation 
of the lithosphere and asthcr.osphcre following an earthquake. Near the fault there is a postscismic 
stress rise whose magnitude can be a significant fraction of the stress reduction due to the earth- 
quake. This stress recovery is due to a change in the rigidity of the lithosphere and is accompanied 
by only small changes in strain. By contrast the stress changes associated with asthenospheric flow 
mimic the strain changes. The amplitude of these changes is small compared to the peak values of 
the stress changes associated with lithospheric relaxation, but the spatial extent of the changes is 
broader. 

The paper has also considered the dependence of coseismic and postscismic surface displace- 
ments, strains, and stresses on fault depth, height, and length (strike-slip faults) and on lithospheric 
thickness. Among the findings are tiiat the peak values of the deformations and stresses increase 
with increasing fault height and decrease with fault depth. Distant deformations and stresses in- 
crease with fault area but are relatively insensitive to the source depth so long as the depth 
remains sufficiently shallow. Values of the postscismic deformations and stresses associated with 
asthenospheric flow are increased by educed values of lithospheric thickness. 
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FIGURE CAPTIONS 


Figure 1 , Viscoelastic Model for Deformation and Stress Due to an Earthquake Along a 
Strike-Slip Fault 

Figure 2. Shear Stress, <r 12 (t). Fault extends along x axis from -L to +L. Tick mark 
spacing = 2L. 

Figure 3. Normal Stress, a u (t) 

Figure 4. Normal Stress, a 2 z(t) 

Figure 5a, Shear Stress, a l2 (t). Versus Distance From Fault, y 

5b. Shear Strain, e 12 (t), Versus Distar.ee From Fault, y 
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Figure 1. Viscoelastic Model for Deformation 


and Stress Due to an Earthquake Along a Stnke-Slip Fault 




Figure 2. Shear Stress, o ]2 (t). 


Fault extends along x axis from -L to +L Tick mark ; 











